Effects of Smooth Boundaries on Topological Edge Modes in Optical Lattices 
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Since the experimental realization of synthetic gauge fields for neutral atoms, the simulation of 
topologically non-trivial phases of matter with ultracold atoms has become a major focus of cold 
atom experiments. However, several obvious differences exist between cold atom and solid state 
systems, for instance the small size of the atomic cloud and the smooth confining potential. In 
this article we show that sharp boundaries are not required to realize quantum Hall or quantum 
spin Hall physics in optical lattices and, on the contrary, that edge states which belong to a smooth 
confinement exhibit additional interesting properties, such as spatially resolved splitting and merging 
of bulk bands and the emergence of robust auxiliary states in bulk gaps to preserve the topological 
quantum numbers. In addition, we numerically validate that these states are robust against disorder. 
Finally, we analyze possible detection methods, with a focus on Bragg spectroscopy, to demonstrate 
that the edge states can be detected and that Bragg spectroscopy can reveal how topological edge 
states are connected to the different bulk bands. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices provide a unique 
experimental setup for studying properties of solid state 
systems in a very clean and well controlled fashion (JJ [2] . 
Particularly interesting in this context is the experimen- 
tal implementation of artificial gauge fields for neutral 
atoms [TJ [3]-[5][6], simulating for instance time-reversal 
symmetry breaking magnetic fields [7HTT] [T2] or a cou- 
pling of the atom's internal spin degree to its angular mo- 
mentum [T3UT6] . The realization of these effects will open 
a path for precise simulations of a large class of topolog- 
ically non-trivial systems such as quantum Hall (QH) or 
quantum spin Hall (QSH) phases. Creating topological 
states of matter with cold atoms is particularly attrac- 
tive because of the precise control of physical parameters 
such as the hopping amplitude and interaction strength, 
allowing the possibility to observe strongly interacting 
topological phases in lattice experiments. However the 
implementation of artificial gauge fields for neutral atoms 
is only one experimental challenge in simulating topolog- 
ical phases in optical lattices [T7U2Q] . Experiments must 
overcome the difficulties provided by the finite size of the 
lattice and the soft boundary of the system, caused by 
a trapping potential that is smoothly varying in space. 
Finite size leads to a finite overlap of spatially separated 
counterpropagating edge states and therefore to possible 
backscattering processes, decreasing the robustness of the 
edge states against external perturbations [21j[22]. While 
this is not a very serious restriction for optical lattice po- 
tentials, which are relatively pure, the effects of the soft 
boundary of the optical lattice system may significantly 
change the properties of the edge states characterizing 
topological insulators in finite systems. Whereas recent 
publications identify the soft boundaries as an unwanted 
restriction or propose how to avoid them by implement- 
ing artificial sharp boundaries to their system [23], we 
demonstrate in this article that soft boundaries will lead 
to interesting additional features, either not present or at 



least not visible in systems with sharp boundaries. For 
this purpose, we investigate different trap shapes and ge- 
ometries, which are realizable in optical lattices and dis- 
cuss their specific influence on the cold atom system. 

This article is organized in the following way. First, in 
Sec. [ll| we present the theoretical model under consider- 
ation, a QH Hamiltonian in the tight-binding approxima- 
tion for spin polarized fermions confined in an additional 
trapping potential. In Sec. |III[ we present our results for 
the stripe geometry, discussing in detail the properties 
of the edge states in systems with a hard wall boundary, 
a harmonic trap and a quartic trapping potential. In 
Sec. [TVl we study the shape of the edge states in a com- 
pletely trapped system and investigate the suitability of 
several detection methods as tools to probe the system 
experimentally, including Bragg spectroscopy. Finally, in 
Sec. [V] we provide some conclusions. 

II. THE MODEL 

The model we consider is similar to the ones proposed 
in [7J [11] , experimentally realizing time-reversal symme- 
try breaking topological insulators with ultracold atomic 
gases. This model describes a two-dimensional (2D) sys- 
tem of spin-polarized fermionic atoms subjected to a 
square optical lattice, experiencing an artificial Abelian 
gauge field A that induces an artificial uniform magnetic 
field perpendicular to the lattice, B = Be z , which is 
similar to the celebrated Hofstadter Hamiltonian [24] on 
the lattice. In our system, the gauge field A enters the 
first-quantized Hamiltonian of the system in form of the 
minimal coupling p — >• p — ^A, which leads to the Hamil- 
tonian: 

H = (p - -A) 2 /2m + W(x) + U(x). (1) 

c 

Hamiltonian ([I]) contains the optical lattice potential 
W and a spatially dependent scalar potential V which 
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allows for the inhomogeneity of the lattice, caused by 
the finite width of the laser beams creating the lattice 
or additional external potentials such as a harmonic trap 
or an artificial hard wall boundary. For the moment we 
leave the detailed shape of V arbitrary, and only assume 
that the non-local matrix elements of V are negligible 
(i.e. (Z|V|ra) = 5j jm (Z|V|Z), where (l\ is the Wannier state 
at lattice site Z), which is reasonable in our case since the 
potential is either varying slowly compared to the lattice 
spacing a, or is a step function. 

The second quantized form of Hamiltonian ([I]) in the 
tight binding approximation then reads 

H = -tJ2 4e i2 ^'-c m + J2 Vic\ Cl . (2) 

l,m I 

The operator c\ here denotes the fermionic creation op- 
erator at lattice site Z, with its respective annihilation 
operator c t . The first term is the well known nearest 
neighbor (NN) hopping with amplitude t, and is com- 
plex due to the Peierls phases 27r0/ ?m that are a re- 
sult of the gauge field [23] . The second term corre- 
sponds to the inhomogeneity V with the local matrix 
elements V\ = (Z|V|Z). The phases <^ m = -j^ J*™ A • d\ 
are not uniquely defined by the magnetic field and de- 
pend on the gauge chosen. In this paper, we choose the 
common Landau gauge A = (0, £?£,()), which leads to 
<t>l i7n = a ■ xi(S yij y m +i - S yil y m -i), where x t and y t are the 
coordinates of lattice site I with lattice spacing a = 1 and 
a = represents the flux per plaquette in units of the 
magnetic flux quantum, $o = h/e. Setting e = h = 1, 
we obtain a = ^- for the square lattice. Throughout the 
rest of this article we choose the hopping t as the natural 
energy unit of our system. 

In the following sections we will restrict our analysis 
to the case where a = 1/6 or a = 2/5 respectively. Our 
results for these two cases can easily be generalized to 
other cases where a = p/q, with q G N, and where 
topological edge states are predicted [25 . 

The experimental realization of a similar model was 
proposed in [23], where the authors consider a spinful 
fermionic system subjected to an artificial gauge field 
that simulates a magnetic field of the form B = Bae z , 
where a = ±1 is the spin quantum number. This model 
preserves time-reversal invariance and therefore allows 
for the realization of QSH phases in optical lattices. Be- 
cause of the time-reversal symmetry, our analysis also 
applies to this model when spin remains a good quan- 
tum number, and we will mention the corresponding 
QSH phases throughout the text. So far, we have not 
accounted for a Zeeman splitting due to an external mag- 
netic field, a spin-orbit coupling or a staggering poten- 
tial, all realizable in optical lattices [25] . The physics 
caused by these additional effects are indeed very inter- 
esting and leaving them out may seem quite restrictive, 
but the results we discuss in this article are quite general 
and require only that the states are topological and do 
not rely on the detailed nature of the edge states. 



III. EDGE STATES IN CYLINDRICAL 
GEOMETRIES 

The defining property of topological insulators in a 
semi-infinite system is the emergence of gapless edge 
states which are localized at one edge and robust against 
perturbations of the system, e.g. potential or magnetic 
disorder. Furthermore, the presence of these states is 
the origin of the currents measured in QH [25-27] and 
QSH samples [28-30 which are well known to be strictly 
quantized when S z is a good quantum number. Topolog- 
ical phases are typically distinguished by the transport 
properties of the edge states, specifically by the quantized 
charge (or mass for neutral atoms) that is transported at 
a single edge [5TJ [32] . One method of determining the 
topological quantum number for a given system is there- 
fore to calculate the energy spectrum on a cylindrical 
geometry and to evaluate the transport properties of the 
edge states directly. Alternatively, one may determine 
the topological quantum numbers from the dispersion re- 
lation of filled bands [35] or the corresponding eigenstates 
[52U34] in the corresponding infinite system. 

In this section we will focus on a cylindrical geometry 
and determine the spectral functions of the system of in- 
terest via exact diagonalization of the Hamiltonian for a 
finite system of size 100 x 100 [36 . We discuss the prop- 
erties of the edge states of the system by analyzing the 
integrated spectral function in quasi- momentum space 
and real space for several kinds of boundaries and show 
the robustness of the edge states against perturbations 
by switching on a disordered potential. 

A. Identification of topological invariants 

Topological phases can be characterized either by an- 
alyzing the band structure properties of the infinite sys- 
tem, or in terms of the transport properties of the system 
in a confined geometry. While the first approach is insen- 
sitive to the specific shape of the confining potential, the 
latter may in principle strongly depend on these details. 
In this section we discuss the properties of edge states at 
an infinite wall boundary, realized by open boundary con- 
ditions at the edges of the cylinder, henceforth referred to 
as stripe geometry. For this kind of boundary the topo- 
logical quantum number of the infinite system is equiv- 
alent to the transport coefficient v of the finite system, 
a relation known as the bulk-boundary correspondence 
[34] [35] . The coefficient v counts the difference in number 
of forwards-moving and backwards-moving states at the 
Fermi edge, which represents the net transport for low- 
energy excitations and hence the quantized edge current 
Ie [25]. Explicitly we have [37] 

v™ = ^2sign(dk y e am (k y )), (3) 

where a m labels the states at the Fermi edge with en- 
ergy e am (ky) = €f and m = L,R for the left and 
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FIG. 1. (Color online) Integrated spectral function for a system described by Q with flux a — 1/6 and stripe geometry with 
100 x 100 lattice sites [36]. Three experimentally relevant confinements of the form V(x) — Vo(x/L) s are shown: a) hard 
wall, S — » oo, b) quartic confinement, S = 4 and c) harmonic confinement, 6 = 2. The spectra in the upper row, pz,(£; y ,cj), 
show the k y -dependence along the periodic direction, integrated over the left half of the confinement direction, see text, and 
the real space spectra of the lower row, p(x,oj), show the x-dependence along the confinement direction. To the right of the 
figure are the transport coefficients, calculated using |3| with the Fermi edge set to the corresponding dotted line. For hard 
wall and quartic confinement there is an appreciable number of bulk bands and the edge states are clearly distinguishable, 
whereas within harmonic confinement we consider almost all of the states to be edge states. In each case we indicate left edge 
states in red (gray), while the remaining bulk states are shown in black, corresponding to the left half of the cylinder in real 
space, see lower plots. Two possible approaches exist for designation of edge and bulk states in softly confined systems. As 
seen in the upper row, energy regions with well defined topological quantum numbers can be identified in the spectrum. The 
corresponding states can be designated as part of the edge, and the remaining ones as the bulk. Alternatively, we can define 
the edge as the point at which no states have energies within the range of energies covered by states at the center of the trap. 
We use the latter designation, although there is little difference between the two methods. 



right edge, respectively. Eq. ([3| can be obtained by 
applying the well-known Laughlin argument to a cylin- 
drical geometry and subsequently following the proce- 
dure described in [34], where no details of the trapping 
potential are required. For the gauge A = (0,£fo,0), 
the single particle Hamiltonian obeys the symme- 
try if(x, p) = H{— x, — p), which leads to vl = —vr. 
Throughout this article we will only consider the Hall 
transport coefficient for the left edge v = z/£ 5 which 
is identical to the topological Z quantum number of 
the infinite system and determines the Hall conductance 
&xy — ve 2 /h. The topological Z2 quantum number v<i, 
which indicates QSH phases in the corresponding spin- 
1/2 system [23] can then be obtained, if S z is a good 
quantum number, by [38j [39] 

v 2 = \v\ mod 2. (4) 



If V2 = 1 the system will exhibit a QSH phase. 

In this paper we make an explicit distinction between 
an "edge state" and an "edge mode". An edge state al- 
ways refer to an eigenstate of the Hamiltonian that is 
localized to one edge, whereas an edge mode refers to 
a series of edge states that are smoothly connected in 
momentum-space. Although this distinction is not nec- 
essary for hard-wall systems, it is required for soft bound- 
ary systems. 



B. Cylindrical geometry with open boundary 
conditions 

We first consider a system described by Q with a step 
potential V that is zero for \x\ < L x /2 and infinite else- 
where, with L x sufficiently large, providing a hard- wall 
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FIG. 2. (Color online) Integrated spectral function pL(k y ,uj) 
for a system described by |2} with flux a = 2/5 and hard- wall 
(left) and quartic (right) confinement. Edge states are shown 
in red (gray), while bulk states are shown in black. 



boundary at the edges of the cylinder. 

Since the quasi- momentum in the ^/-direction is a 
well-defined quantum number and we are interested 
in transport coefficients for this direction, it is con- 
venient to represent the spectrum of the system in 
terms of the integrated spectral density pL{k yi uo) = 

f-L/2 dxp(x-, x, k y , cj), where the spectral function is de- 
fined as 



p(x, x r , k y ,Lj) = — 2Im(x, k y 



1 



uo - H + iO+ 



W,k y ), (5) 



We integrate only over the left half of the system in 
real space, so as to separate the left from the right edge 
states [40]. In Fig. [IJi), upper, the integrated spectral 
density pL{k yi uo) is shown for a = 1/6. One can identify 
the bulk states which are grouped into six thick bands 
and the edges states, which close the gaps between the 
bands. To determine the transport coefficients and possi- 
ble topological phases, we place the Fermi edge in a bulk 
gap and apply (|3| to the dispersion of the edge states. 

There are several phases visible in this system. If the 
Fermi edge lies within a bulk band, the system is in a 
trivial metallic phase. If the Fermi edge lies between the 
third and fourth band, there are a set of Dirac points with 
a linear dispersion and the phase is a semi-metal. In the 
bulk gaps, which lie between the other bands, the edge 
modes place the system in a quantum Hall phase with 
v = — 1 and —2 for the first and second bulk gap, respec- 
tively, and inverted for the third, fourth bulk gap. In 
shorthand, we can specify the phases between the bands 
by gap 1 / 6 = { — 1, — 2, D, 2, 1}, where D represents Dirac 
points. In the analogous spin- 1/2 system, v = ±1 in- 
dicates a QSH phase whereas v = ±2 corresponds to a 
normal insulator, due to lack of topological protection. 

In addition we also investigate the case where a = 2/5, 
see Fig. [2 for which we find gap 2 / 5 = {2, —1, 1, —2}. 
Note that there are no Dirac points for a = 2/5. The 
differences between gap!/ 6 and gap 2 /5 appear as differ- 
ent real space behaviors within soft confining potentials 
that are not visible within a hard- wall confinement, as 
we demonstrate in the next section. 



C. Cylinder with soft boundaries 

With soft boundaries, it becomes relevant to look at 
the spectra of the system in real space along the x 
axis, p{x,uj) = j dk y p(x,x,k y ,u), as well as the par- 
tially integrated spectra in quasi-momentum space along 
the k y axis, pL(k yi uj). The quasi-momentum spectra 
allows us to extract transport coefficients and discuss 
the dispersion of the edge modes. On the other hand, 
the real space spectrum exhibits some unusual features 
for different trapping conditions. We consider a lattice 
of size 100 x 100 and trapping geometries of the form 
V(x) = V (x/L) 6 where V = lOt and L = 50a is cho- 
sen such that V(x = 50a) = V(x = —50a) = lOt which 
is larger than the energy spanned by the infinite system 
~ St. Three particular values of 5 are relevant to ex- 
periment: S — ^ 00 which reproduces hard- wall boundary 
conditions, S = 4 for quartic confinement and S = 2 for 
harmonic confinement. In most optical lattice experi- 
ments, the confining potential is a result of the Gaussian 
envelope of the finite beam width of the lasers and in the 
center of the trap we may approximate this confinement 
by its leading order harmonic term. However, it is has 
been suggested [41] that one can remove the harmonic 
term by superimposing an anti-trapping Gaussian beam 
of different detuning to the trapping beam, which then 
promotes the quartic term to the leading order approx- 
imation of the trapping potential, i.e. V(x) oc x 4 . This 
scheme was realized in optical lattice experiments to im- 
prove quantum phase diffusion experiments [42] 03 . We 
investigate these trapping geometries below in further de- 
tail. 



1. General features and preservation of topological 
invariants 

In Fig.[l] we show a comparison between pL(k y ,uj) and 
p(x,uj) for a = 1/6 and hard wall, quartic and harmonic 
confinements that are relevant to experiment. One can 
see that the potential does not gap the system, and edge 
modes continue to connect the bands. The transport 
coefficients of the soft boundary systems, indicated to the 
right of Fig. [I] are insensitive to the trapping potential. 
In other words, there exist energy ranges in which we can 
identify a transport coefficient, which is identical for all 
confinements we consider. 

Comparing the /^-dependent soft-boundary spectra of 
the upper row in Fig. [I] we make two observations: 1) 
we can readily identify highly degenerate regions of bulk 
bands in the hard-wall and quartic confinements, and 
2) we find that the dispersion of edge modes that are 
present within quartic confinement do not change no- 
ticeably when changing the confinement to the harmonic 
trap. In contrast, the rest of the spectrum is significantly 
modified, such that the ratio of bulk to edge states is very 
small. To define such a bulk region in the soft boundary 
system, we assume the edge begins at a distance from the 
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FIG. 3. (Color online) False color diagram of the integrated spectral density pL(k y ,uj) for a system described by |2| with flux 
a = 1/6, stripe geometry, and a confining potential, V(x) = Vo(x/L) s . Bulk bands are indicated in black and edge modes as 
colored curves, also marked as (a), (6), (c), (d). There exist several true crossings and avoided crossings in the spectra which 
combine to preserve the topological invariants for any confinement exponent 5. Auxiliary states of the corresponding edge 
modes are shown with dashed curves. The auxiliary states do not influence the topological phases of the system, since they 
always come in pairs with opposing velocities. 



trap center where none of the states at this point overlap 
in energy with any of the states in the very center of the 
trap (i.e. at x = 0). 

Analogously, we can clearly identify a bulk region from 
the x-dependent spectra in the bottom row of Fig. [I] for 
the quart ic trap, but not in the case of the harmonic trap. 
From this we conclude that the quartic trap is likely the 
best trapping potential for observing effects of both the 
bulk system and topological edge states in an experi- 
mental setup, if it is not feasible to artificially implement 
hard wall boundaries as proposed, for example, in [23] . 
Furthermore, we observe no overlap between states of 
different edges, which has been proposed to destroy edge 
states via couplings between the edges [21]. This again 
shows that the edge states are topologically protected 
and robust against external changes in the potential. 

In Fig. |2l we also show a comparison of pL(k y ^uj) for 
the case where a = 2/5 between a hard wall (left column) 
and a quartic confined system (right column). We again 
calculate the transport coefficients of both systems and 
list these next to the plots to show that these also coincide 
for all trapping potentials. 

To better understand the details of the rather compli- 
cated spectra of the quartic trap, we choose to follow the 
edge modes of the hard-wall confinement by smoothly 
varying the trapping exponent 5. We show plots for 
5 = 60, 16, 8 and a = 1/6 in Fig. [3j where we have ar- 
tificially colored the spectra to indicate each edge mode. 
5 = 60 represents a very steep trap, and is almost identi- 
cal to the hard-wall case: with the color designation, one 
sees that the blue (marked as (b)) and yellow (c) edge 
states are only present in the 2nd and 3rd bulk gaps, 
respectively, whereas the red (a) and green edge states 
(d) span two bulk gaps. As the confinement is made 
softer, we see that an edge mode may cross the BZ more 
than once, and that the energy range of the edge states 
changes, e.g. with 5 = 8 the red state (a) now extends 
into the 3rd bulk gap. However, whenever this occurs, 



the state forms an avoided crossing at some higher energy 
with a different edge state and is forced downwards in 
energy, a process which preserves the value of the topo- 
logical invariants. We represent this in the false color 
diagram by a dotted line for parts of the edge states that 
are non-topological, i.e. not connecting different bulk 
bands. For 5 = 16, 8, we can consistently see this oc- 
curring in the most energetic edge mode (colored green, 
(d)), which extends above the highest bulk band, and 
forms an avoided crossing with the non-topological edge 
state created by the effect of the trapping on the highest 
bulk band. 

Note that, due to the trapping potential, several edge 
states that belong to the same edge mode may exist for 
one value of energy. 

2. Merging and splitting of edge states 

When the number of edge modes changes, as the Fermi 
edge crosses a bulk band in the hard-wall boundary sys- 
tem, either an edge mode must be created, or an edge 
mode must merge into either the bulk band itself or with 
another edge mode. In the soft-boundary system we can 
see some very non-trivial behavior that shows the com- 
plexity of these processes. 

We first focus on the real space spectra of the a = 2/5 
flux system under quartic confinement, see Fig. [4] In 
the lowest gap, we see that two different edge modes, 
which evolve between the first and second bulk band, 
merge into a single edge mode, which evolves between the 
second and third bulk band. In the hard wall system, this 
mode is localized to a single site in the x-direction and 
can only be observed in quasi- momentum space. In the 
quartic trap, the edge states leaving the first bulk band 
follow the shape of the quartic potential and one may 
expect the same for the states leaving the second bulk 
band. As one sees in Fig. [4] this is not the case. The 



states leaving the second bulk band immediately start to 
merge with the edge states from the first bulk band and 
the result is only a single mode at each edge, evolving 
between the second and third bulk band. Although it 
is not possible to determine topological invariants from 
real-space spectra, we can link this merging behavior to 
the ky-sp&ce spectra of Fig. [2] and see that it leads to the 
correct topological quantum number v = — 1. The same 
effect is again observable between the fourth and fifth 
bulk band. Interestingly, the merging of these modes 
does not take place via a simple overlap of the states, 
but a gap in real space with negligible spectral weight 
exists between the states originating from the bands and 
the newly-formed edge mode. 

In the a = 1/6 flux system, we also see the opposite 
effect: the splitting of a single bulk band, to connect edge 
modes of different bands which are energetically well sep- 
arated. In Fig. [I] the integrated spectral density p(x, uS) 
shows that the modes leaving the second and fifth bulk 
bands each split into two curves, where a single eigenstate 
has large amplitudes on two spatially separated lattice 
sites. We interpret this splitting as a process that facil- 
itates the connection between different bands which we 
observe in Fig. [3j For example, the outer part of the 
mode leaving the second bulk band can be seen to merge 
at higher energies with the mode that is a product of 
the third and fourth bands. This connection between 
the bands is analogous to the avoided crossings that we 
observe in the /^-dependent spectra in Fig. [3] This very 
non-trivial behavior of modes within the outer region of 
the system, combined with transport coefficients which 
are identical to the topological quantum numbers, given 
by the transport coefficients of the infinite system, indi- 
cates that the soft edge states are of topological origin. 



To further verify this, we address in section HID the ro- 
bustness of these states against perturbations in terms of 
a disordered background potential. 



3. Relation of edge states and bulk bands 

When we look more closely at the dispersion of the 
edge modes, we can see an interesting connection to 
the bulk bands of the system. We focus on the quasi- 
momentum spectra for the case a = 1/6 shown in Fig. [3] 
for increasing confinement exponent S. 

The dispersion of an edge mode leaving a given bulk 
band can be described on two different quasi- momentum 
scales. For a small range of k y , the dispersion mimics 
that of its associated bulk band, and this behavior be- 
comes more prominent for smaller 5. This can be seen 
for the lowest edge modes, colored red (a) and blue (6), 
e.g. the red (a) mode has a locally flat dispersion, mir- 
roring the flatness of the lowest band. However, when 
avoided crossings have occurred, such as for the yellow 
(c) and green (d) edge modes, the dispersion of an edge 
mode cannot simply be described by one band alone and 
corresponds to a mixture of bands. 
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FIG. 4. (Color online) Integrated spectral density p(x,uj) for 
a system described by Eq. [2]on a 100 x 100 lattice with quartic 
confinement V. In the left figure, the flux is chosen to be a = 
2/5, while in the right figure a = 1/6. In the left figure, the 
edge states leaving the second bulk band immediately merge 
with the edge states of the first bulk band to form a state 
spatially localized between the bulk bands with the correct 
transport coefficient. In the right figure, the states leaving 
the second bulk band split up, i.e. they localize to more than 
one point in space. The inner part of these states merges with 
the edge states of the third band at higher energies. Similar 
behavior is observed for the 3rd, 4th and 5th band. As pointed 
out in the text, this non-trivial behavior is an indication of 
the topological origin of the edge states. 



On the other hand, considering the Brillouin zone as 
a whole, the edge modes become more flat in momen- 
tum space, the smoother the confining potential is in 
real space. This flattening is a direct result of the num- 
ber of accessible sites at the edge. The number of lattice 
sites, n e d ge , that are available for an edge state between 
e.g. the first and the second bulk band, is the number 
of sites i that fulfill e\ — eo J$ V{xi) ;$ 62 — eo, where 
ei (62) are the maximum (minimum) energies of the first 
and second bulk band, respectively, and eo is the mini- 
mum energy of the first bulk band [22 . In the hard wall 
system, n e d ge = 1 but becomes larger the smoother the 
confining potential becomes in real space. An interest- 
ing result for the a = 1/6 flux per plaquette, is that the 
flatter the potential becomes, the flatter the lowest gap 
edge modes become, with a corresponding increase of the 
effective mass of system's excitations: 



m 



( d 2 



00. 



(6) 



This is generally true for soft confinements, as pointed 
out in [22], but the edge state structure in the case a = 
1/6 allows this feature even for relatively steep potentials. 



D. Robustness of soft edge states in stripe 
geometries 

One of the most important properties of topological 
edge states is their robustness against even large pertur- 
bations, which leads to clearly detectable quantized Hall 
conductance in impure experimental setups. In optical 
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FIG. 5. (Color online) Integrated partial spectral density 
pL(k x ,oj) of a system described by Eq. [2]on a 120 x 60 lattice 
with a — 1/6. Right: with an additional binary disordered 
potential, given by Eq. [7] and A max = 0.5t. Left: with the 
disordered potential being set to zero. The robust edge states 
are still clearly pronounced and gapless, while the former bulk 
bands are smeared out and show a mobility gap (not shown 
here but obtainable from the Anderson-localized bulk eigen- 
states) . 



lattice experiments which are, by construction, very clean 
realizations of condensed matter Hamiltonians, perturba- 
tions such as disorder are usually not an issue. However, 
since disordered potentials can be implemented in a con- 
trolled manner [44-46 , it is of interest to thoroughly in- 
vestigate how robust the edge states are against these 
kinds of perturbation. Here, we address this question 
for soft boundaries. The general argument, which illus- 
trates the robustness of edge states in condensed matter 
systems, is the lack of possible backscattering processes 
[32] . Counter-propagating edge states are localized on 
opposite edges of the system and are very well separated 
spatially. Therefore, in huge condensed matter systems, 
these states have vanishing spatial overlap and backscat- 
tering from impurities is completely suppressed. In con- 
trast, in finite systems different edge states from opposite 
edges will have a finite overlap in real space, which the- 
oretically allows for backscattering processes, and there- 
fore disorder may lead to the opening of a gap in the 
spectrum. However, as we will see from our numerical 
results, even in very small systems (~ 60 lattice sites in 
x-direction) this effect is not observable. 

To verify the robustness of the soft edge states nu- 
merically, we perturb the system described by Eq. [2] by 
adding a disordered background potential 



v dis , 



i4 c i 



(7) 



where A^ is distributed randomly, either by a binary dis- 
tribution, G {0,A max }, V7, or by a uniform distribu- 
tion, Aj G [0, A max ], V7 [47 . For all realizations, we found 
that the edge states stay robust and still connect the dif- 
ferent bulk band regions without opening a gap up to 
disorder strengths of about A max w 0.5£ for binary disor- 
der and even larger strengths for uniform disorder. 

For example, in Fig. [5] the integrated partial spectral 
density pL{k y ,uS) is shown for a 120 x 60 lattice system 



with uniformly distributed A/ and A max = 0.5t. There 
is clearly no gap in the spectrum and although quasi- 
momentum is no longer a good quantum number, the 
edge states in momentum space are very sharply centered 
around a particular value of k y and remain delocalized in 
the ^/-direction as they were for the system without disor- 
der [48]. In contrast, some of the bulk states now consist 
of many quasi-momentum components (not shown in our 
figure) and therefore become localized to a region much 
smaller than the system size, which can be termed An- 
derson localization. 

We have also addressed larger systems with larger 
boundary regions. These systems contain more and more 
edge states in a given bulk gap, which may possibly lead 
to different backscattering processes between edge states 
located at the same edge and therefore open gaps in the 
spectrum after disorder is introduced. To exclude these 
possibilities, we studied system sizes of up to 60 x 240 
lattice sites without finding any indication of gaps in the 
spectrum or localization of the edge states up to disorder 
strengths of A max = OM. 



IV. DETECTION METHODS 

So far we focused on a semi-infinite system with stripe 
geometry. However, realistic systems in optical lattice ex- 
periments are confined to a finite region in all dimensions 
by the finite beam width of the lasers. In this section, 
we determine the spatial wave functions of a 2D system 
trapped in both the x- and ^/-directions and discuss pos- 
sible detection methods of the resulting edge states. 



A. Eigenstates of the completely trapped system 

We determine the eigenstates of a system with a con- 
fining potential V that varies in the x- and ^-directions 



V(x,y) = V 



(z)'+(2)' 



(8) 



The parameter S determines the shape of the trap and 
the possible eigenstates. For 6 — > oo the system is again 
confined by hard walls in both directions, while for 5 = 4 
and 5 = 2 the system is in a quartic and harmonic con- 
finement, respectively. For harmonic confinement, we ex- 
pect the eigenstates that are extended over various lattice 
sites to be circularly symmetric, whereas in the quartic 
case the potential is no longer circularly symmetric and 
the states take on a shape that is sometimes referred to 
as a squircle [49]. For this analysis we again will restrict 
ourselves to a system with 100 x 100 lattice sites and 
a trapping potential with a minimum value of Vq = lOt 
along each edge of the lattice. We again focus on a = 1/6. 

In Fig. [6] we clearly see the different real space distri- 
bution of edge states compared to bulk states. The bulk 
states are delocalized over a region of about 30 x 30 lat- 
tice sites, while the edge states for a given energy follow 



FIG. 6. (Color online) Wave function \^(x,y)\ 2 as a function 
of the lattice spacing a of different eigenstates of the system 
within complete quartic confinement. The real space spectral 
density shown to the left is a cross section of the complete 
system: p(x, y = 0, lj). Three particular eigenstates have been 
shown, and their energies indicated by arrows to the spectral 
density. The wave functions belong to a) a bulk region and b) 
and c) to a pair of edge states splitting up after leaving the 
second bulk band. 

the isolines of the quartic potential, and are strongly con- 
fined to these regions. Comparing figures b) and c) for 
the completely trapped system one can see the splitting 
of the edge states leaving the second bulk band as the two 
states have a weak overlap with one another. The shape 
of the edge states in the quartic confinement looks similar 
to that which one would expect in the hard wall system 
(like those explicitly shown in [50]) and differs only at 
the very corners of the system. We therefore expect sim- 
ilar single particle excitations for the quartic confinement 
as for the hard wall confinement when probing the edge 
states in experiment. 

The situation slightly changes when looking at the har- 
monically confined system. There, the confining poten- 
tial is circularly symmetric and one may expect that the 
eigenstates reflect this symmetry. The wave functions of 
the harmonically confined system are shown in Fig. [7J As 
already seen from the spectral density plotted in Fig. [T] 
the former bulk region is tightly confined to very few lat- 
tice sites in the center of the trap, which makes it diffi- 
cult to define a bulk region in the harmonic trap. On the 
other hand, the edge states chosen reflect the radial sym- 
metry of the trapping potential and again are localized 
along the isolines of the trapping potential. This already 
indicates that for the harmonic confinement we expect 
very different excitation dynamics than for the quartic 
and hard wall confinement, where significant parts of the 
eigenstates are quasi-one-dimensional. 

For harmonic confinement, one can solve the contin- 
uum model analytically in the absence of the lattice 
[5TJ [52] and the resulting wave functions are quite similar 
to those from the lattice calculation. The major differ- 
ence in the continuum case is that no edge states from 
different bulk bands merge, since the hall conductivity 



FIG. 7. (Color online) As in figure [6] but with complete har- 
monic confinement. The three states are shown from a) a bulk 
region and b) and c) edge states of states belonging to differ- 
ent bulk bands. Spatial coordinates in units of the lattice 
spacing a. 

always increases by one when passing a bulk band and 
the different bulk bands are not connected. Additionally, 
angular momentum is only a good quantum number in 
the continuum case. 



B. Bragg spectroscopy 

An important question concerning topological non- 
trivial phases in ultracold atoms is whether the edge 
states are detectable with existing experimental tools. 
Due to the lack of stationary transport in optical lattice 
experiments, it is not feasible to directly measure the 
Hall conductance, and one has to consider alternative ap- 
proaches [22j [53] . Several possibilities for detecting edge 
states or topological quantum numbers in optical lattice 
experiments have been proposed. Some require careful 
experimental implementation such as Bragg [54] [55] or 
Raman spectroscopy [56 and others take advantage of 
easily accessible observables like time-of-night (ToF) pat- 
terns or density profiles. 

Density profile measurements were proposed by Umu- 
calilar et al. [57J to directly separate the bulk and edge 
densities between different bands. However, as already 
pointed out in [22], these profiles do not show the re- 
quired structure, as can be seen in Fig. [I] (lower): the 
bulk bands all occupy approximately the same real space 
extent. Hence, this method is not applicable to topolog- 
ical systems in general. 

Alternatively, ToF measurements have been proposed 
by Zhao et al. [58] to exhibit minima and maxima that 
depend on the topological number of the system. While 
this is true for the specific cases they were investigating 
and also for our system in the case of a = 1/6, we found 
that it is not valid in the case of a = 2/5 and therefore 
cannot be reliably used as a detection method in exper- 
iment. In contrast, Alba et al. [59] propose using ToF 
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measurements as a method to identify skyrmions, by fo- 
cusing on topological properties of pseudo-spin vectors 
within the Hamiltonian on the Bloch sphere. However, 
this method focuses on bulk properties rather than the 
edge modes that we consider here. 

We choose to focus instead on Bragg spectroscopy, 
which probes the dynamical structure factor S(q,uj) of 
the underlying system. Bragg scattering of topological 
insulators in optical lattices has been previously con- 
sidered for the case of the quantum anomalous Hall 
effect [54]. However, no inhomogeneity of the lattice 
was considered. Recently, Goldman et al. [55 have in- 
vestigated Bragg spectroscopy theoretically, considering 
shaped lasers to probe angular momentum states within 
circularly symmetric traps. While this is a novel imple- 
mentation to enhance the detection of edge states, we 
demonstrate that one is able to observe edge states using 
a simple linear Bragg coupling which, due to technical 
limitations, may be the only option available to a par- 
ticular experiment. Furthermore one can observe differ- 
ences in Bragg spectroscopy between the various bands 
that we show is not due to chirality considerations. We 
do not propose an explicit experimental setup and sim- 
ply assume that one can measure the dynamical struc- 
ture factor directly. One such proposal to measure this 
precisely in an optical lattice is the so called "shelving 
method" [55] . 

When performing Bragg spectroscopy, the system is 
illuminated by two laser beams, described by wave vec- 
tors Pi,P2 and frequencies uj\ = pic, uj^ = respec- 
tively, and the differences in these quantities, q = P1—P2 
and = c^i — c^2, allows for transitions between differ- 
ent eigenstates of the original system. The Hamiltonian 
describing the interaction of the system with the laser 
beams is then given by 



e-^ t * t (q + p)*(p) + h.c.) 



(9) 



where is a field operator, creating a particle with 

real momentum p and T is the coupling strength of the 
lasers [6QH62] , 

The dynamical structure factor in linear response the- 
ory for an infinite homogeneous system is directly con- 
nected to the density-density correlation function Xq,q(u) 

S(q,uj) = --lm Xq , q (uj) (10) 

7T 

via the fluctuation dissipation theorem [63]. For our case, 
we have to evaluate S(q,uj) for the inhomogeneous sys- 
tem, where the quasi-momentum is no longer a good 
quantum number. Within the linear response approxi- 
mation, and accounting for the finite size of the system 
and finite time of the measurement process, we find: 



S(q,co) = \T\ 2 Aj2 



n x {l - n^Ax^q) | 5 



/Li, A 



■w A ) 2 + A 2 ' 



(11) 



respectively, of the state v. We introduce a Lorentzian 
broadening factor A, to allow evaluation in a system of 
finite size. The scattering amplitude A\ jfJL (q) is the prob- 
ability of a particle in state \i to scatter into the state A 
by gaining momentum q and is given by the integral 



A 



d 3 re-^C(r)^A(r). 



(12) 



Here, A, \i label the single-particle eigenstates of the sys- 
tem, n v and uj v are the occupation number and energy, 



After determining the single particle eigenstates of the 
system, we can directly calculate the dynamical struc- 
ture factor. Because we are focusing on the detection of 
edge states, we investigate a system with a Fermi energy 
located in a bulk gap at ep = —2t (see Fig. [I]), where 
an edge state is located. There are now four general 
scattering processes possible, edge^edge, edge^bulk, 
bulk^edge and bulk^bulk. Edge^edge scattering is 
clearly distinguishable by analyzing the dynamical struc- 
ture factor. Given a frequency lj, the set of possible mo- 
mentum transfers allowed to another edge state is very 
limited because the edge states are well localized in mo- 
mentum space. For the case of edge^bulk scattering, 
many different momenta are accessible and therefore we 
see a signal regardless of the value of q. This means for 
a fixed momentum transfer qr, S(q,uj) as a function of 
uj consists of a £-peak approximately around uj = qvp 
[64] and a smeared out region, where the bulk bands are 
located. This can be seen in Fig. [8] (left), where the first 
peak indicates edge— ^edge scattering and the second and 
third peaks correspond to edge— )>bulk and bulk^bulk 
scattering. 

On the other hand, for a fixed frequency cj, the re- 
sponse in momentum space describing edge— ^edge scat- 
tering looks quite different from that obtained from 
edge-^bulk scattering, as one can see from Fig. [8] a) 
and b), respectively. For the quartic confinement, the 
edge states form squircles in real space (Fig. [6|, which 
means that low energy excitations are most favorable in 
x- or ^/-direction, resulting in the square-like structure of 
S(q,uj) in Fig. [8j which is approximately described by 
{^allowed} = {(q x ,q y ) I max{|^|, \q y \} « q = v f /uj). In 
contrast, the dynamical structure factor of bulk^bulk 
scattering from the first to the second band is smeared 
out and depends on the Fermi surface of the occupied 
level at ep and the Fermi surface of the unoccupied band 
€p = €p + uj. The allowed momenta are approximately 
described by {<? a i lowed } = {(q x ,q y ) I \q x \ + \q y \ ~ € F /v F } 
and form a rough square which is rotated by tp = 7r/4 
compared to the edge— )>edge scattering. Note that as a 
result of the structure of A^\(q), where a minimal spatial 
overlap of the two spatial wave functions is needed for 
obtaining a finite scattering amplitude, high frequency 
edge^edge scattering is exponentially suppressed be- 
cause the presence of the trap causes energetically sep- 
arated edge states to be localized to different distances 
from the center of the trap. This does not occur in the 
equivalent hard-wall system. 

For the harmonically confined system, S(q,uj) as a 
function of uj for fixed q is qualitatively the same as 
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FIG. 8. (Color online) Dynamical structure factor for a sys- 
tem with quartic confinement and Fermi energy in the first 
bulk gap (cf = —2t). Left: S(q,uj) for fixed momentum q 
as a function of uj. The first peak belongs to edge— ?>edge 
scattering and its position is sensitive to q and can be writ- 
ten, for small uj, as uu q — VF,edg e q- The second and third 
peak belong to edge— )>bulk scattering from edge states into 
the third and fourth bulk bands, located around e = and 
to bulk— )>bulk scattering from the first to second bulk bands, 
where the frequency is independent of q. No signal appears 
of scattering from edge states to the second bulk band, lo- 
cated at uj = 0.5t, indicating a disconnection between these 
states, i.e. these states have vanishing matrix elements of the 
Bragg operator. Right: S(q,uj) for fixed frequency as a func- 
tion of momentum transfer q for a) edge— ^edge scattering at 
uj = 0.2£, b) bulk— )>bulk scattering processes at uj = 1.5t. 



in the quartic system. In contrast, S(q,uu) as a func- 
tion of q for fixed uj for edge— )>edge scattering looks 
quite different than for quartic confinement. As seen 
in Fig. [7| the edge states have a circular symmetry and 
therefore no momentum transfer direction is preferred, 
which leads to the set of allowed states forming a cir- 
cle {^allowed} = {(qx,Qy) I 

Fig. H 



q 2 y 



j/vp}, shown in 



An important discovery of our calculations is that there 
is an obvious absence of spectral weight at frequencies 
where we expect signals of edge^bulk scattering. To 
highlight this, we calculate an artificial Bragg response 
where we allow only initial states in the energy range 
— 2.5£ < e < — 2t for transitions to states of higher energy. 
This means any signal due to possible bulk^edge or 
bulk^bulk transitions is suppressed. The spectra shown 



in Fig. [10] demonstrate edge^edge signal for the first 
bulk gap, edge— ^bulk signal to the third and fourth bulk 
bands, but conspicuously absent signal for the edge— )>2nd 
bulk band transitions, which would be expected for Bragg 
frequencies 0.5 < uj < 1.5. This implies that the first 
edge and second edge/bulk are disconnected, i.e. have 
a vanishing matrix element of the Bragg operator. It is 
possible to predict this behavior from the dispersion of 
the edge states (see Fig.|3|, as one can see that the lowest 
edge mode, colored in red (a), passes unimpeded through 
the second bulk band, and never displays an avoided 
crossing with the blue (b) edge mode of the second band 




(units of 7r/a) -1-1 



/ a ) % -1 -1 (units of tt/o) 

(units of ir/a) 



FIG. 9. (Color online) S(q,uj) for a fixed frequency uj = 0.2t 
as a function of momentum transfer q, for a Fermi energy 
€f = — 2t. Left: S(q,uj) for the quartic confined system. 
The system shows a strong response when one component of 
q — (qx,q y ) has an absolute value \q x , y \ = qo = oj/vf because 
excitations along the x-axis, y-axis are most favorable (see 
Fig. [6]). Right: S(q,uj) for the harmonically confined system. 
Here, the response is close to circularly symmetric in g-space, 
reflecting the shape of the eigenstates. No particular direction 
is anymore favorable, as long as the absolute value of \q\ = qo 
is fixed. 
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FIG. 10. (Color online) Dynamical structure factor S(q,uj) 
for fixed momentum q as a function of frequency, for a Fermi 
energy cf = — 2t, as seen in Fig. [8] but with artificially 
suppressed bulk— >bulk scattering processes. The first peak 
belongs to edge to edge scattering processes and is sensi- 
tive to the momentum transfer q with approximate frequency 
uj q = VF,ed ge q- The broadened peaks around uj = 1.5t belong 
to edge— )>bulk scattering to the third and fourth bulk bands. 
It is clearly visible that there is no scattering from the edge 
states to the second bulk band, which is located at uj = 0.5t. 



or the second band itself, while it always merges with the 
3rd or 4th band (with which we find non-vanishing ma- 
trix elements of the Bragg operator). For higher energies, 
and strong confinements, we see the opposite behavior of 
avoided crossings between red (a) edge modes and yellow 
(c) edge modes, indicating that one can expect a finite 
Bragg response from transitions between these states. 
Note that the lack of edge^bulk scattering is not a 
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result of the soft-boundaries inhibiting real-space over- 
lap. We have performed equivalent hard- wall boundary 
calculations where real-space overlap is guaranteed but 
we again observe an absence of signal for disconnected 
edge— >bulk transitions. Note also, that we do not expect 
to observe a clear signal for large frequency edge^edge 
transitions regardless of the type of trap, as there is larger 
range of states beneath the Fermi-edge that can be ac- 
cessed with the Bragg laser. Hence, many different values 
of k y will contribute, leading to a blurred signal. 

V. CONCLUSION 

In this article, we analyzed the properties of 2D topo- 
logical edge states in softly confined systems with a con- 
finement in one direction of the form V(x) = Vo(x/L) s . 
By varying the confining potential from a hard-wall to a 
quartic or harmonic potential, we showed that the topo- 
logical properties of the edge states in specific bulk gaps 
do not depend on the steepness of the confining potential, 
while a confinement sharper than harmonic is required to 
achieve an appreciable bulk region of the lattice. We sug- 
gest that quartic confinement is suitable to observe both 
edge state and bulk properties, which may be realized by 
superimposing attractive and repulsive Gaussian beams. 
Furthermore, we observed the emergence of robust aux- 
iliary edge states, which provide additional structure to 
edge modes but do not influence the topological quan- 
tum numbers. The main feature of these auxiliary states 
is that they connect edge states which are spatially sep- 
arated to bulk bands of the system. This provides a 
mechanism to preserve topological invariance, as soon as 
the edge states and bulk bands become spatially sepa- 



rated. In these cases the band-structure exhibits a series 
of avoided crossings that act to preserve the topological 
invariant. An analysis of the spectral density of softly 
confined systems in real space revealed the splitting and 
merging of edge states from different bulk bands, which 
is also indicative of their topological nature. 

We also determined the wave functions of eigenstates 
in a completely trapped system and showed how these 
depend on the confining potential. With these, we cal- 
culated the dynamical structure factor which can, for in- 
stance, be measured by Bragg spectroscopy. We found 
that the dynamical structure factor can reveal the edge 
and bulk states of the system and their overlap. 

In summary, we demonstrated that topological 
properties in ultracold atomic systems with artificial 
gauge fields are not sensitive to the trapping potentials 
available in optical lattice experiments and that the edge 
states of these systems can be clearly detected via Bragg 
spectroscopy. We believe that soft boundaries provide 
more detailed insight into the behavior of edge states, 
which cannot be observed in hard- wall systems, and are 
therefore worth investigation in their own right. 
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